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1. Introduction 

Why yet another review of big bang nucleosynthesis (BBN)? Given that 
the basic physics of the synthesis of the light elements D, 3 He, 4 He and 7 Li 
in the big bang was thoroughly discussed over 25 years ago ||, ^ it is nat- 
ural to wonder why so many papers on the subject continue to appear on 
a regular basis. The reason for this is two-fold. First, as observational de- 
terminations of light clement abundances have advanced, the situation has 
become more instead of less uncertain. In particular it has become apparent 
that inferring the primordial values of the abundances from contemporary 
observations is fraught with uncertainty, given our limited understanding 
of the chemical evolution of galaxies. Fortunately observers have risen to 
the challenge and developed sophisticated techniques to look further back 
into our past, at nearly pristine primordial material. However large dis- 
crepancies, most likely of a systematic nature, have subsequently emerged 
between different determinations, so the improvement in precision has not 
led to an increase in accuracy! Nevertheless this is a healthy development 
in that it has provided a refreshing perspective on the strong claims made 
in the past decade concerning the consistency of BBN predictions with the 
inferred primordial abundances, and on the stringent constraints thus in- 
ferred on new physics. The second development concerns recent efforts to 
improve the accuracy of the theoretical predictions of the abundances, and, 
just as important, to quantify the uncertainties. In this talk I will discuss 
these issues and summarize their implications for the dark matter problem 
and for new physics, viz. N v ^ 3. Other recent reviews [|[ ||, || may be 
consulted for a different perspective than mine || 
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2. Theoretical Calculations 



We now have excellent semi-analytic insights into the physics of BBN which 
enable the 4 He abundance to be calculated correctly to within ~ 1 — 2% Q 
and the 'left-over' abundances of D, 3 He, and 7 Li to be estimated to within 
a factor of ~ 2 — 3 ||. However for precision work, and in particular, to 
determine the uncertainties, it is necessary to use the Wagoner computer 
code H which has been improved and updated with the latest values of the 
nuclear cross-sections and made freely available to the community [[[o| . 

There have been two important developments on this front. First the 
rates for the weak interactions which determine the ncutron-to-proton ra- 
tio have been carefully calculated beyond the Born approximation, with 
allowance for all zero and finite temperature radiative, Coulomb and finite 
(nucleon) mass corrections [[l2[ fL3| . With further allowance for finite tem- 
perature QED corrections to the equation of state of the plasma and the 
residual coupling of the neutrinos to the plasma during e + e~ annihilation, 
the 4 He abundance has been computed with an accuracy of ~ 0.4% jl2| , 
where the dominant source of uncertainty comes from the present experi- 



mental determination of the neutron lifetime: r„ = 885.3 ±2.0 s 14 . This 
is an impressive feat and it is believed that all relevant physical processes 
have now been consistently taken into account. 

The uncertainties in the abundances of the other light elements are of 
course considerably greater and are, moreover, correlated with each other. 
The standard practice has been to use the Monte Carlo method to sample 
the error distributions of the relevant reaction cross-sections which are then 
input into the numerical code, thus enabling well-defined confidence levels 
to be attached to the theoretically predicted abundances (l5| [ll| . However 
this is computationally expensive and, moreover, needs to be repeated each 
time any of the input parameters are changed or updated^] To overcome 
this handicap we have developed a method based on linear error propaga- 
tion which requires the numerical code to be run just once to determine the 
covariance (or error) matrix; simple % 2 statistics can then be used (rather 
than maximum likelihood methods as with Monte Carlo simulations) to 
determine e.g. the best-fit value of r\ = n^/n^, the nucleon-to-photon ra- 
tio fTi| . This method has recently been extended to consider departures 
from the standard model, viz. an effective number of neutrinos N u during 
BBN different from 3, so that correlated limits on rj and N v can be ex- 
tracted for any given set of input abundances |Q. The results agree well, 
where comparison is possible, with similar exercises using the Monte Carlo 
plus maximum likelihood method |o). All the calculations have been 
encoded in a simple code which is available from a website jll| . 

1 For example, it has recently been argued that uncertainties in the cross-sections 
for some key nuclear reactions are in fact smaller than were estimated earlier lla] . 
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The method is briefly as follows. First define the four relevant elemental 
abundances as Y2 = D/H, Y3 = 3 Hc/H, Yj = 7 Li/H, and finally Y4 as the 
mass fraction of 4 He (often termed Yp). These depend both on the model 
parameters 77 and N v , and on a network of nuclear reactions Rk'. 

Y i = Y i (r,,N v , {R k }) . (1) 

The rates for the 12 essential nuclear reactions which have to be considered 
plus their associated uncertainties have been discussed in detail Eg] a nd we 
use these, apart from the updated value of the neutron lifetime |H . Now 
for a small change of the input rate Rk (Rk —* Rk+S Rk), the corresponding 
deviation of the i-th elemental abundance (Y — > Yi + 5Yi) as given by linear 
propagation reads 

5Y i ( V ) = Y i (r 1 )^2x ik (r ] ) S -^ , (2) 

k 

where the functions \k (rj) represent the logarithmic derivatives of Yj with 
respect to R k : 

XM = m^RW)- (3) 

In general, the deviations SYi are correlated, since they all originate from 
the same set of reaction rate shifts {8 Rk}. The global information is con- 
tained in the error matrix (also called covariance matrix) p2] ], which is a 
generalization of the "error vector" 5Yi in eq. (|^). In particular, the abun- 
dance error matrix af^r]) obtained by linearly propagating the input ±ler 
reaction rate uncertainties ±A Rk to the output abundances Y reads: 

oj j (v)= Y i(v)Y j (v)J2 X M X M (^) 2 ■ (4) 

This matrix completely defines the abundance uncertainties. In particular, 
the lcr abundance errors ai of Yj are given by the square roots of the 
diagonal elements, 

ViW = V^SW . (5) 



while the error correlations pij can be derived from eqs. (|J,|5|) through the 
standard definition 

PijW = —ri — rr ■ 6 

We have checked that the propagation of errors is indeed linear and then 
computed the error matrix above, finding good agreement with the results 
obtained by Monte Carlo [pL ITrl . 
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Next we must input the inferred values of the primordial abundances, 
Yi±Ui. Taking the errors Ui in the determinations of different abundances 
Y i to be uncorrelated, the experimental squared error matrix afj is simply 

w ij = SijViVj , (?) 

where Oy is Kronecker's delta. The total (experimental + theoretical) error 
matrix Sfj is then obtained by summing the matrices in eqs. (^,0): 

S?-(»7)=4(r?)+4 • (8) 
Its inverse defines the weight matrix Wij(rj): 

WiM = [s$ j (vT l - ( 9 ) 

The x 2 statistic associated with the difference between theoretical ( Yj) and 
observational (Yl) light element abundance determinations is then |22j| : 

x 2 fa) = E K fa) ~ ^1 ' w v fa) ' K fa) - ^ ■ (io) 

ij 

Contours of equal % 2 can then be used to set bounds on the parameters 
r\ and N v at selected confidence levels. For standard BBN we set = 3 
and minimization of the \ 2 then gives the most probable value of 77, while 
the intervals defined by \ 2 = Xmin + ^ X 2 gi ve the likely ranges of 77 at the 
confidence level set by A \ 2 ■ 

Below we have provided a table of polynomial fits to the central val- 
ues of the abundances, Yj = ao + a±x + a^x 1 + a^x 3 + 04a; 4 + 0,5a; 5 , with 
x = log 10 (?7/10~ 10 ) in the range 0-1 jL7|. The value of Y4 from the Wag- 
oner code was corrected as described earlier [ 6| and is in excellent agreement 
with a subsequent precision calculation Jl2 |. These fits, along with simi- 
lar fits to the logarithmic derivatives Xik(x) | p7[ , have been encoded in a 
javascript program which plots the abundances with associated uncertain- 
ties and allows the corresponding values of r] to be read off for a given input 
abundance |E3[. We now proceed to discuss the observational situation. 
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3. Observed Abundances 



The abundances of the light elements synthesized in the big bang have been 
subsequently modified through chemical evolution of the astrophysical en- 
vironments where they are measured p4| . The observational strategy then 
is to identify sites which have undergone as little chemical processing as 
possible and rely on empirical methods to infer the primordial abundance. 
For example, measurements of deuterium (D) can now be made in quasar 
absorption line systems (QAS) at high red shift; if there is a "ceiling" to 
the abundance in different QAS then it can be assumed to be the pri- 
mordial value, since it is believed that there are no astrophysical sources 
of D. The helium ( 4 He) abundance is measured in H II regions in blue 
compact galaxies (BCGs) which have undergone very little star formation; 
its primordial value is inferred either by using the associated nitrogen or 
oxygen abundance to track the stellar production of helium, or by simply 
averaging over the most metal-poor objects J25fl , (We do not consider 3 He 
which can undergo both creation and destruction in stars |24j and is thus 
unreliable for use as a cosmological probe.) Closer to home, the observed 
uniform abundance of lithium ( 7 Li) in the hottest and most metal-poor 
Pop II stars in our Galaxy is believed to reflect its primordial value p6| . 

However as remarked earlier, improvements in observational methods 
have made the situation more, instead of less, uncertain. Large discrepan- 
cies, of a systematic nature, have emerged between different observers who 
report, e.g., relatively 'high' [Eft [2£| or 'low' ||(J values of deuterium in 
different QAS, and 'low' @ §§ or 'high' §| |§ values of helium in BCG, 
using different data reduction methods. It has been argued that the Pop II 
lithium abundance p5[ |36| , |37| may in fact have been significantly depleted 
down from its primordial value |3^, Q , with some observers arguing to the 
contrary 

Rather than take sides in this matter we instead consider several combi- 
nations of observational determinations, which cover a wide range of possi- 
bilities, in order to demonstrate our method and obtain illustrative best-fits 
for j] and N v . 



3.1. Deuterium 

The first observations of D absorption due to a QAS at redshift z = 3.32 
towards Q0014+813 suggested a relatively 'high' value of J27]. ps| 

Y 2 = 1.9 ±0.4 x 10" 4 , (11) 

and was consistent with limits set in other QAS. However subsequent in- 
dependent measurements in QAS at z = 3.572 towards Q1937-1009 and at 
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z = 2.504 towards Q1009+2956 found a much lower abundance of |3(| 



Y 2 = 3.4 ±0.3 x KT 5 . (12) 

More recently, observations of a QAS at z = 0.701 towards Q1718+4807 
have also yielded a high abundance of ps{ 

% = 3.3 ± 1.2 x 10~ 4 . (13) 



These discrepant measurements may result from spatial variations of 
the D abundance due to localized astrophysical processes |p8| or due to 
inhomogeneous nucleosynthesis [fl2) . More prosaic explanations invoke sys- 
tematic effects. For example, new data on the Q0014+813 system J43| 
reveal the presence of an 'interloper' hydrogen line contaminating the deu- 
terium absorption feature. It has also been suggested Q that the discor- 
dance between the 'high' and 'low' values may be considerably reduced if 
the analysis of the H+D profiles accounts for the correlated velocity field 
of bulk motion, i.e. mesoturbulence, rather than being based on multi- 
component microturbulent models. It is then found that a true primordial 
abundance of 

T 2 = (3.5-5.2) x 10~ 5 (14) 

is compatible simultaneously (at 95% C.L.) with the 'low' and 'high' values 
quoted above in eqs. ( p^ , [T^ ) if the data on these objects is appropriately 
analysed. Whether this is indeed the correct explanation in all cases is 
however not clear. Obviously more measurements are needed and at the 
moment all the suggested values of the primordial deuterium abundance 
quoted above merit consideration. 



3.2. Helium 

The primordial helium abundance is found to be ]32] | 

%= 0.234 ±0.002 ±0.005 , (15) 

from linear regression to zero metallicity in a set of 62 BCGs, based largely 
on observations which gave a relatively 'low' value J3l[ . However using data 
on a sample of 27 BCGs from an idependent survey, a significantly higher 
value of 

Ti = 0.243 ± 0.003 , (16) 

has been derived from a new analysis which uses the helium lines them- 
selves to self-consistently determine the physical conditions in the H II 
region Q and specifically excludes those regions (with anomalously low 
He I line intensities) which are believed to be affected by underlying stellar 
absorption. In particular it has been demonstrated that there is strong 
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underlying stellar absorption in the NW component of I Zw 18, which was 



included in earlier analyses 31 . Excluding this, the average of the values 



found in the two most metal-poor BCGs, I Zw 18 and SBS 0335-052, is |J 

Yl = 0.245 ±0.004 . (17) 

Again further work is needed to resolve the situation but we will consider 
both the 'low' and 'high' values for the primordial helium abundance. 

3.3. Lithium 

The controversy here concerns whether the observed approximately con- 
stant 'Spite plateau' in the abundance of 7 Li in Pop II stars in the halo 
reflects its true primordial value or has been depleted down from an ini- 
tially higher abundance. There is agreement that the average value in the 
hottest and most metal-poor stars is about [^5|, 

Yr = 1-6 ±0.36 x 10" 10 , (18) 

but there is dispute about whether there is significant scatter suggestive of 
some depletion mechanism, or a slight trend of decreasing 7 Li abundance 
with increasing metallicity suggestive of secondary cosmic ray production. 
A recent analysis based on 41 stars finds no evidence of any depletion or 
post-BBN creation and suggests a slightly higher primordial abundance 

Y 7 = 1.73 ±0.21 x 10~ 10 . (19) 

(The systematic error may have to be raised by ~ 50% to allow for the 
uncertainty in the oscillator strengths of the lithium lines |26| .) 

However there are hot, metal-poor Pop II stars which appear identical 
to those which define the 'Spite plateau' but which have no detectable 7 Li 
|j5| , [§8| ], indicating that some depletion mechanism does exist. Recently 
it has been argued that the observed abundance has been depleted down 
from a primordial value of |4j| 

%= 3.9 ±0.85 x 10" 10 , (20) 

the lower end of the range being set by the presence of the highly overde- 
pleted halo stars and consistency with the 7 Li abundance in the Sun and 
in open clusters, while the upper end of the range is set by the observed 
dispersion about the 'Spite plateau' and the 6 Li/ 7 Li depletion ratio. A 
somewhat smaller depletion is suggested by other workers who find a 
primordial abundance of Y-j = 2.3 ± 0.5 x 10 -10 . 

The implications of all these possibilities need study. 
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4. Implications for r\ and baryonic dark matter 
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Elsewhere |18| we have discussed various combinations of the above data 
sets and provided a handy computer programme jllj to enable the reader 
to consider other permutations/combinations. Here we will present results 
for just two possible (although mutually incompatible) selections of mea- 
surements which we name data set "A" and data set "B" : 



(21) 



(22) 

(In subsequent figures we also consider a slight variation of data set B in 
which we use eq. ( |l7| ) instead of eq. ( [l6| ) for the 4 He abundance.) 

Figure |l| shows the theoretical predictions compared with the data. We 
calculate and minimize the x 2 an d obtain the 95% C.L. ranges allowed by 
each data set by cutting the curves at A% 2 — X 2 ^ Xmin = 3.84: 

Data set A : r) = 1.78±[^ 4 x 10" 10 , (23) 
Data set B : 77 = 5.13±^6 x 10_1 °- ( 24 ) 

The range of 77 for data set A agrees very well with the 95% C.L. range 
obtained independently with the same input data but using the Monte 
Carlo plus maximum likelihood method fll9| . However it is inconsistent 
with a further constraint on 77 coming from a recent analysis of the Lya- 
"forest" absorption lines in quasar spectra. The observed mean opacity of 
the lines requires some minimum amount of neutral hydrogen in the high 
redshift intergalactic medium, given a lower bound to the flux of ionizing 
radiation. Taking the latter from a conservative estimate of the integrated 
UV background due to quasars yields the constraint |l6| 

77 > 3.4 x lO" 10 . (25) 

Using Q^h 2 ~ 77/2. 73 x 10~ 8 this corresponds to S1n > 0.0125/7 -2 , where 
f?N is the density of nucleons in ratio to the critical density and h is the 
present Hubble parameter in units of 100 km s -1 Mpc -1 . 

Thus data set B is favoured by the Lya-forest constraint. The implied 
nucleon density is significantly higher than the nucleonic matter in stars 
and X-ray emitting gas in groups and clusters today Q , 

fi N = — ~ 0.0025/1- 1 + 0.0046/1- 1 - 5 , (26) 

Pc 

implying that most of the nucleons are presently in some dark form. 



Standard BBN predictions vs data (1a errors) 




x = 109,0(77 X10 10 ) 

Y 2 = D/H 

Y 4 = 4 He mass fraction 
Y 7 = 7 Li/H 

Data set A 

Y 2 = (1.9±0.4)X10~* 
Y 4 = 0.234±0.0054 
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Figure 1. Standard BBN predictions (dotted lines) in the planes defined 
by the abundances Yz, Y4, and Yr fFflj . The theoretical uncertainties are 
depicted as lcr error ellipses at x = log 10 (?//10" 10 = 0, 0.1, 0.2, . . . , 1. The 
crosses indicate the two observational data sets (with la errors). 



5. New Physics and the Possibility of N„ 7^ 3 

We should also take into account that the goodness of fit can be affected 
by a non-standard Hubble expansion rate during BBN, e.g. due to the 
presence of new neutrinos. Although the Standard Model contains only 
N v = 3 weakly interacting massless neutrinos, the recent experimental 
evidence for neutrino oscillations Q may require it to be extended to 
include new superweakly interacting massless (or very light) particles such 
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as singlet neutrinos or Majorons. These do not couple to the Z° vector 
boson and are therefore unconstrained by the precision studies of Z° decays 
which establish the number of SU{2)\ J doublet neutrino species to be [jl4j 

N„ = 2.993 ± 0.011 . (27) 

However, as was emphasized some time ago ]49[| , such particles would boost 
the relativistic energy density, hence the expansion rate, during BBN, thus 
increasing the yield of 4 He. This argument was quantified for new types 
of neutrinos and new superweakly interacting particles [ ]50| in terms of 
a bound on the equivalent number of massless neutrinos present during 
nucleosynthesis : 



^ 2 \T L 

% 

where gi is the number of (interacting) helicity states, /b = 8/7 (bosons) 
and /f = 1 (fermions), and the ratio TijT v depends on the thermal history 
of the particle under consideration ]51|] , For example, Ti/T v < 0.465 for 
a particle which decouples above the electroweak scale such as a singlet 
Majoron or a sterile neutrino. However the situation may be more compli- 
cated, e.g. if the sterile neutrino has large mixing with a left-handed dou- 
blet species, it can be brought into equilibrium through (matter-enhanced) 
oscillations in the early universe, making Ti/Tv ~ 1 ]5^ |. Moreover such 
oscillations can generate an asymmetry between v e and v Cl thus directly af- 
fecting neutron-proton interconversions and the resultant yield of 4 He ||^] . 
This can be quantified in terms of the effective value of N v parametrizing 
the expansion rate during BBN, which may well be below 3! Similarly, 
non-trivial changes in N v can be induced by the decays p4[ or annihila- 
tions J55| of massive neutrinos (into e.g. Majorons), so it is clear that it is 
a sensitive probe of new physics. 

The values of N v and rj are correlated since the effect of a faster ex- 
pansion rate can be compensated for by the effect of a smaller nucleon 
density. Indeed increasingly stringent lower bounds on r\ inferred from 
assumed upper limits to the primordial deuterium abundance (which were 
deduced from chemical evolution arguments) have been used to set increas- 
ingly stringent upper bounds on N v . These have ranged from 4 downwards 
|3(| , culminating in one below 3 which precipitated the so-called "crisis" for 
standard BBN [f57| , and was interpreted as requiring new physics. How- 
ever as cautioned before |5S|| , there are large systematic uncertainties in 
such constraints on N v which are sensitive to our limited understanding 
of galactic chemical evolution. Moreover it has been emphasized that 
the procedure used earlier |5(| to bound N v was statistically inconsistent 
since, e.g., correlations between the different elemental abundances were 
not taken into account. This can of course be addressed by the Monte 



10 



Carlo method, using which it was shown J39[ that the conservative ob- 
servational limits on the primordial abundances of D, 4 He and 7 Li allow 
N v < 4.53 (95% C.L.), significantly less restrictive than earlier estimates. 
However because the Monte Carlo method is computationally expensive 
and the exercise needs to be redone e.g. whenever the input abundances 
change, it would clearly be advantageous to extend our linear error propa- 
gation approach to include N„ as a parameter. 

This is, in principle, straightforward, since it simply requires recalcula- 
tion of the functions Yi and \k at the chosen value of N v . Since it would 
be impractical to use extensive tables of polynomial coefficients, we have 
devised some formulae which, to good accuracy, relate the calculations for 
arbitrary values of N v to the standard case N v — 3, thus reducing the 
numerical task dramatically [Q. 

As is known from previous work j^], the synthesized elemental abun- 
dances D/H, 3 He/H, and 7 Li/H (i.e., Y%, Y3, and Yj in our notation) are 
given to a good approximation by the quasi-fixed points of the correspond- 
ing rate equations, which formally read 

dY 

^oci^yxyxHT, (29) 
+,- 

where the sum runs over the relevant source (+) and sink (— ) terms, and 
{<jv)t is the thermally-averaged reaction cross section. Since the temper- 
ature of the universe evolves as dT /At oc —T 3 ^/g^, with the number of 
relativistic degrees of freedom, g+ = 2 + (7/4)(4/ll) 4 / 3 ^ (following e+e" 
annihilation) , the above equation can be rewritten as 

l^-i^E^^HT, (30) 

which shows that Y 2 , Y3, and Y7 depend on ij and N v essentially through 

1 /2 

the combination r\j gj . Thus the calculated abundances Y2, Y3, and Yj (as 
well as their logarithmic derivatives Aj&) should be approximately constant 
for ^ 

log rj — — log 5* = constant , (31) 

as we have verified numerically 

This suggests that the values of Yi and of \k for N v = 3 + AiV^ can be 
related to the case N v = 3 through an appropriate shift in x: 

Y l {x,Z + AN U ) ~ Yi(x + aAN v ,3) , (32) 
A lfc (x,3 + AN,) ~ \ tk {x + Cl AN„,3) , (33) 



where the coefficient a is estimated to be ~ —0.03 from eq. ( |3l| ) (at least 
for small AN„). In order to obtain a satisfactory accuracy in the whole 
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range (x, N v ) € [0,1] x [1,5], we in fact allow upto a second-order vari- 
ation in AN V , and for a rescaling factor of the Y^s. As regards the 4 He 
abundance, a semi-analytical approximation Q also suggests a relation be- 
tween x and N v similar to eq. (|3l|), although with different coefficients, and 
indeed functional relations of the kind above work well also in this case. 
However, in order to achieve higher accuracy and, in particular, to match 
exactly the result of the recent precision calculation of I4 which includes 
all finite temperature and finite density corrections (T^[ , we also allow for a 
rescaling factor for the A^'s. All these calculations have been embedded in 
a compact Fortran programme available from a website JTT[ , which allows 
easy extraction of joint fits to x and N v for a given set of elemental abun- 
dances, without having to run the full BBN code, and with no significant 
loss in accuracy. We now proceed to discuss these fits. 



6. Joint fits to rj and N v 

Figure |^ shows the abundances for various N v compared for illustration 
with data set A. While there is consistency between theory and data for 
x ~ 0.2 - 0.4 and N„ = 3, for N u = 2 {N v = 4) the Y 2 data prefer values 
of x lower (higher) than the Y4 data. Therefore, we expect that a global 
fit will favor values of (x,N u ) close to (0.3,3). Similarly, a fit to data set 
B should favour values of (x,N v ) ~ (0.7, 3) as seen in figure ||. 

To quantify this we show in figure || the results of joint fits using the 
abundances of data set A. The abundances Y 2 , Y i: and Y 7 are used sep- 
arately (upper panels), in combinations of two (middle panels), and all 
together, without and with the Lya-forest constraint^] on 77 (lower panels). 
In this way the relative weight of each piece of data in the global fit can 
be understood at glance. The three C.L. curves (solid, thick solid, and 
dashed) are defined by x 2 — Xmin = 2-3, 6.2, and 11.8, respectively, corre- 
sponding to 68.3%, 95.4%, and 99.7% C.L. for two degrees of freedom (77 
and N v ), i.e., to the probability intervals often designated as 1, 2, and 3 
standard deviation limits. The x 2 is minimized for each combination of 
Yi, but the actual value of Xmin ( an d the best-fit point) is shown only for 
the relevant global combination Y% + Y4 + Y7(+Lya). For this data set, 
the helium and deuterium abundances dominate the fit, as can be seen by 
comparing the combinations Y 2 + Y 4 and Y2 + Y4 + Y7. The preferred values 
of x are relatively low, and the preferred values of N v range between 2 and 
4. Although the fit is excellent, the low value of x is in conflict with the 
Lya-forest constraint on 77, as indicated by the increase of Xmin from 0.02 
to 8.89. 

2 This bound is not well-defined statistically but we can parametrize it for example 
through a penalty function Xi^yaiv) = 2-7 X (3.4 X 1Q~ 10 /r/) 2 which can be added to the 
X 2 derived from our fit to the elemental abundances. 
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Data set 
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Figure 2. Primordial abundances I4 ( 4 He mass fraction) Y2 (D/H) and 
Y'j ( 7 Li/H), for N u = 2, 3, and 4. Solid and dashed curves represent the 
theoretical central values and the ±2<r bands, respectively. The grey areas 
represent the ±2a experimental bands for the data set A M . 

As shown in figure ^| there is no such problem with data set B which 
favors high values of x because of the 'low' deuterium abundance. The 
combination of Y2 + Y7 isolates, at high x, a narrow strip which depends 
mildly on N„. The inclusion of Y4 selects the central part of such strip, 
corresponding to N v between 2 and 4. As in figure ||, the combination Y4 + 
Y'j does not appear to be very constraining. The overall fit to Y2 + Y4 + Y 7 
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Data set 



Allowed bands 



y 2 = (3.4 ± 0.3) x io- 5 :::::::::::::: theory (± 2a) 

Y 4 = 0.245 ± 0.004 

Y 7 = (1.73 ± 0.21) x IO"' DATA (± 2a) 

N, = 2 N„ = 3 H„ = 4 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



log(r/io) iog(r/, ) l°g(?7io) 



Figure 3. Same as in figure |, but for data set B Q. 

is acceptable but not particularly good, mainly because I2 and Yj are only 
marginally compatible at high x. On the other hand, the Yz + 3*4 + Yj 
bounds are quite consistent with the Lya- forest constraint. 

Elsewhere [[l8) we have presented results for other possible combinations 
of input abundances. As shown in figure ^ (cf. (20|), one can also consider 
orthogonal combinations to those above, e.g. 'high' deuterium and 'high' 
helium, or 'low' deuterium and 'low' helium. The latter combination implies 
N v ~ 2, thus creating the so-called "crisis" for standard nucleosynthesis 
. Conversely, the former combination suggests N v ~ 4, which would 
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Data set '. 

Y 2 = (1.9 ± 0.4) x 10"" 
Y 4 = 0.234 ± 0.0054 
Y 7 = (1.6 ± 0.36) x 10"' 



Allowed regions '. 



68.3% C.L. 
95.4% C.L. 
99.7% C.L. 



N, 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



N, 



N, 



Y 2 +Y 4 




Y 4 +Y 7 





0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



Y 2 +Y 4 +Y 7 



0.02 



- 4 

- 3 
2 

J 1 



Lya 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

log(?] 10 ) logfaio) logfoio) 



Figure 4. Joint fits to x = log(r?io) and N v using the abundances of data 
set A [||). The abundances Y 2 , I4, and Yj are used separately (upper 
panels), in combinations of two (middle panels), and all together, without 
and with the Lya-forest constraint on rj (lower panels). 



also constitute evidence for new physics. Allowing for depletion of the 
primordial lithium abundance to its Pop II value, relaxes the upper bound 
on N v further, as was noted earlier [p9[ . 
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Data set '. 

Y 2 = (3.4 ± 0.3) x 10~ 5 

Y 4 = 0.245 ± 0.004 

Y 7 = (1.73 ± 0.21) x 10" 



Allowed regions '. 



68.3% C.L. 
95.4% C.L. 
99.7% C.L. 



N, 



N, 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



N, 



Y 2 +Y 4 +Y 7 



3.71 



- 4 

- 3 
2 

J 1 



Lya 



Y 2 +Y 4 +Y 7 +Lya 



- A - 

- 3 - 
2 - 

J 1 Li. 



4.94 



0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

log (77,0) log (77,0) log(?7io) 



Figure 5. Same as in Fig. Q but for the data set B [H 



7. Conclusions 



It is clear from the above discussion that the present observational data on 
the primordial elemental abundances are not as yet sufficiently stable to 
derive firm bounds on rj and N u . Different and arguably equally acceptable 
choices for the input data sets lead to very different predictions for 77, and to 
relatively loose constraints on N v in the range 2 to 4 at the 95% C.L. Thus 
it may be premature to quote restrictive bounds based on some particular 
combination of the observations, until the discrepancies between different 
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Figure 6. The 68% (solid) and 95% (dotted) likelihood contours for the 
number of neutrino species and the nucleon-to-photon ratio, for all four 
combinations of the 'high' and 'low' deuterium and helium abundance mea- 
surements in data sets A and B. 



estimates are satisfactorily resolved. Our method of analysis provides the 
reader with an easy-to-use technique 1 1 to recalculate the best- fit values 
as the observational situation evolves further. 

In order to improve the theoretical predictions further it is necessary 
to know the nuclear reaction rates better. To determine which reaction is 
largely responsible for the uncertainty in a particular elemental abundance 
(for any given nucleon density) is clearly important in this regard. In our 
approach this can be easily determined by "perturbing" the values of the 
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Figure 7. Individual contributions of different reaction rates Rk to the 
uncertainties in Y2, 14, and I7, normalized to the corresponding total errors 
CT2, <T4, and (77, for 77 = 1.78 x 1CP 10 , the best-fit value for data set A |17y . 
Each arrow corresponds to the shift SYi induced by a +ler shift of Rk- 
Some small error components have not been plotted. 



input reaction rates and observing their effect on the predicted abundances, 
namely one can study the contribution to the total uncertainty <7j of Yi 
induced by a +ler shift of R k : 

R k —> Rk + A R k Y -> Yt + SYi ■ (34) 

We choose two particular values of r/, 1.78 x 1CP 10 and 5.13 x 1CP 10 , corre- 
sponding to the best-fit values for data sets A and B respectively and show 
in figures and || the deviations SYi (normalized to the total error c^) 
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Big Bang Nucleosynthesis — Error Components 

at 7] = 5.13X10" 10 




-1 -0.5 0.5 1 -1 -0.5 0.5 



6T 4 /a 4 (SY 2 /a 2 

Figure 8. Same as in figure]]], but for 77 = 5.13 x 10~ 10 , the best-fit value 
for data set B . 

induced by +la shifts in the i?fc's, plotted in the same set of planes as used 
for figure @. The la error ellipses shown in these figures are obtained by 
combining the deviation vectors SYijui in an uncorrelated manner. Several 
interesting conclusions can be drawn from this exercise. As expected, the 
uncertainty in the weak interaction rate R\ has the greatest impact on I4 
for the high value of 77 (figure ^) , since essentially all neutrons end up being 
bound in 4 He. However at the lower value of 77 (figure the uncertainty 
in i? 2 — the "deuterium bottleneck" — plays an equally important role as 
i?i in determining Y4 because nuclear burning is less complete here than 
at high rj. Similarly with reference to the reaction rates -R7, i?io — R12 
which synthesize 7 Li, at low rj it is the competition between R7 and i?n 
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which largely determines Y7, while at high 77 it is the competition between 
i?io and i?i2- The anticorrelation between Y4 and Y2 is driven mainly by 
i? 2 at low 77 and, to a lesser extent, by R4 and _R 5 , while the reverse is 
the case at high 77. The anticorrelation between Y4 and I7 at low 77 is also 
basically driven by i?2, while the correlation at high 77 is due to both R2 
and i?4. Thus we have a direct visual basis for assessing in what direction 
the output abundances Yi are pulled by possible changes in the input cross 
sections Rk- This should be helpful in determining experimental strategies 
to determine key reaction rates more precisely. 

However one might ask what would happen if these discrepancies re- 
main? We have already noted the importance of an independent constraint 
on 77 (from the Ly a- forest) in discriminating between different options. 
However, given the many assumptions which go into the argument [4q] , 
this constraint is rather uncertain at present. Fortunately it should be 
possible in the near future to independently determine 77 to within ~ 5% 
through measurements of the angular anisotropy of the cosmic microwave 
background (CMB) on small angular scales |pQ|, in particular with data 
from the all-sky surveyors MAP and PLANCK [|6l| . Such observations will 
also provide a precision measure of the relativistic particle content of the 
primordial plasma. Hopefully the primordial abundance of 4 He would have 
stabilized by then, thus providing, in conjunction with the above measure- 
ments, a reliable probe of a wide variety of new physics and astrophysics 
which can affect nucleosynthesis (for recent work see |i2[ |64|, |6^, |6(| |67|] ) . 
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